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Abstract 

By using the exact partition function of su(m\n) Polychronakos spin chain associated 
with Atv-i root system, we study some statistical properties of the related spectrum. It is 
found that the corresponding energy level density satisfies the Gaussian distribution and 
the cumulative distribution of spacing between consecutive energy levels obeys a certain 
'square root of a logarithm' law. 
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Exactly solvable one dimensional quantum integrable spin chains and dynamical mod- 
els with long-range interaction [1-11] have attracted much attention in recent years due 
to their appearance in a wide range of subjects like fractional statistics [TJl [32], quan- 
tum electric transport in mesoscopic systems [14], Yangian quantum group [15-17], SUSY 
Yang-Mills theory and string theory [18-20]. Among quantum spin systems with long- 
range interaction, the well known spin-- Haldane-Shastry (HS) model is introduced in 
an attempt to construct an exact ground state which would coincide with the U — > oo 
limit of Gutwiller's variational wave function for the Hubbard model [HE]. A natural 
su(m) generalization of this exactly solvable HS model is constructed by using the 'spin' 
exchange operator associated with the fundamental representation of su(m) algebra [6], [T] . 
Subsequently, it is realized that such HS spin chain may be reproduced from the trigono- 
metric spin Calogero-Sutherland model by applying the 'freezing trick', which basically 
uses the fact that spin and dynamical degrees of freedom of the latter model decouple 
from each other for large values of the coupling constant. Furthermore, a new quantum 
spin chain with long-range interaction is obtained by applying this freezing trick to the 
case of rational spin Calogero model [8]. Lattice sites of this spin chain, which are in- 
homogeneously distributed on a line, are determined through the zeros of the Hermite 
polynomial [10]. This quantum integrable as well as exactly solvable spin system is usu- 
ally known as Polychronakos or Polychronakos-Frahm spin chain in the literature. Both 
HS and Polychronakos spin chains admit su(m\n) super symmetric extensions, in which 
each site is occupied by either one of the m type of bosonic states or one of the n type 
of fermionic states [21-27]. Such supersymmetric spin chains play a role in describing 
some strongly correlated systems in condensed matter physics, where holes moving in 
the dynamical back ground of spins behave as bosons, and spin- 1/2 electrons behave as 
fermions [281129] . 

It should be noted that, the above mentioned HS and Polychronakos spin chains (along 
with their supersymmetric extensions) are all related to the Aat-i type of root system, for 
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which the interaction between any two spins depends only on the difference of their site 
coordinates. Variants of these spin chains associated with other root systems, endowed 
with more general type of interactions, have also been studied in the literature [30-36]. 

By applying the method of freezing trick, it has become possible to compute the 
exact partition functions of su(m) Polychronakos and su(m) HS spin chains associated 
with A N _i root system [HI EI], su(m\n) supersymmetric extensions of these spin chains 
[23l [25] . and variants of such spin chains related to other root systems [33-36]. These 
exact partition functions have turned out to be a very efficient tool for studying some 
statistical properties of the related energy spectra like level density distribution and dis- 
tribution of spacing between consecutive energy levels. It is found that, for sufficiently 
large number of lattice sites, the energy level density of this type of spin chain follows the 
Gaussian distribution with high degree of accuracy [33-38,25]. It is also observed that, 
distribution of spacing between consecutive energy levels for such integrable spin chains 
is not of Poisson type, as may be expected due to a well-known conjecture of Berry and 
Tabor [39]. However it appears that, even though the exact partition function of su(m\n) 
supersymmetric Polychronakos spin chain associated with A^-i root system has been 
derived by using the freezing trick [23], statistical properties of the related spectrum like 
distributions of energy level density and spacing between consecutive energy levels have 
not been analyzed till now. The purpose of this letter is to study these spectral properties 
of the supersymmetric Polychronakos spin chain by using its exact partition function. 

The Hamiltonian of the su(m\n) supersymmetric Polychronakos spin chain associated 
with the An~i root system is given by [23] 




where the zeros of the iV-th order Hermite polynomial, and the supersymmetric 
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exchange operator -Pj™'™' 1 is defined as 

m+n 

Here creation-annihilation operators like Cj a and C^ a are bosonic for a G {1,2, . . . ,m} 
and fermionic for a G {m + 1, m + 2, . . . , m + n}, and these operators act on a restricted 
Hilbert space satisfying the constraint: 

m+n 



— 1 ) 



a=l 

for all j. It may be noted that for the special case like m ^ 0, n = or m = 0, n ^ 0, 
i.e. when all degrees of freedom are either bosonic or fermionic, 7i(' m l n ) in Eq. (fTJ) reduces 
to the non-supersymmetric Polychronakos spin chain, for which statistical properties of 
the spectrum have been studied very recently [38]. In this article we shall analyze the 
spectrum of 7i( m ' n ) in the purely supersymmetric case, when both m and n take nonzero 
values. 

As shown in Ref. [22], the supersymmetric exchange operator ([2]) can be mapped to 
an 'anyon like' representation of the permutation algebra on an appropriate spin space. 
By using this mapping and subsequently applying the freezing trick, the exact partition 
function of su(m\n) Polychronakos spin chain has been derived as [23] 

z {mln \ q ) = V (g)jv • q^ 1 2 (3) 

N f~? n M (q)aMa 2 ---(q)a m -(q) b 1 (q)b 2 ■■■((!) b n ' 

2^=1 a i + l^j = l b 3= N 

where a^'s and b^s are non-negative integers, q = e~ l l kT and the notation: (g)jy 
; I — g)(l — q 2 ) ■ ■ ■ (1 — q N ) (with (q) = 1) is used. Exploiting a connection between 
this partition function and supersymmetric Schur polynomials, it is found that all energy 
eigenvalues of su(m\n) Polychronakos spin chain can be expressed in terms of 'motifs', 
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which are represented by arrays of binary digits like '0' and '1' [23]. For a spin chain 
with iV number of lattice sites, each motif is represented by an array of N — 1 number 
of binary digits. Thus we may write such a motif as (5i, S 2 , ■ ■ ■ , 5jv-i), where 8{ G [0, 1]. 
These motifs characterize a class of irreducible representations of Y(gl(m\n)) Yangain 
quantum group, which span the Fock space of su(m\ri) supersymmetric Polychronakos 
spin chain. Since Y(gl(m\n)) Yangain is the symmetry algebra of this spin chain, all of 
its eigenstates associated with any particular irreducible representation or corresponding 
motif yield the same energy eigenvalue. One can express the energy eigenvalue of such 
degenerate eigenfunctions associated with the motif S = (Si, 62, ■ ■ ■ , Sn-i) as 

AT-1 

m = £ j *j ■ (4) 

3=1 

It should be noted that, in contrast to the non-supersymmetric case, there is no restriction 
in arranging l's or 0's in a motif associated with the supersymmetric Polychronakos 
spin chain [23]. Consequently, some additional motifs and energy levels appear in the 
spectrum of this supersymmetric spin chain in comparison with its non-supersymmetric 
counterpart. In fact, for this supersymmetric spin chain with N number of lattice sites, 
one can construct 2 N ~ X number of motifs by filling up iV — 1 positions arbitrarily with 
either or 1. By using Eq. (J4j) , it is easy to check that the motif (0, 0, ... , 0) yields the 
ground state with zero energy and the motif (1, 1, . . . , 1) yields the highest excited state 
with energy N(N — l)/2. Moreover, depending on the choice of 5, E{$) in Eq. (j3]) can 
take any integer value within the range to iV(iV — l)/2. Consequently, the spectrum of 
su(m\n) supersymmetric Polychronakos spin chain would be equispaced in nature. 

Since Z^ n \q) in Eq. ([3]) represents the partition function of a finite spin system 
with integer eigenvalues, it should be expressed as a polynomial function of q. Indeed, 
with the help of a symbolic software package like Mathematica, we can explicitly write 
down Z^ n \q) as a polynomial of q for wide range of values of m, n and N. If the 
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term q Ei appears in such a polynomial, then E± will represent an energy eigenvalue and 
the coefficient of q Ei will determine the degeneracy factor corresponding to this energy 
level. Let us denote this degeneracy factor or 'level density' associated with the energy 
level Ei as D^ m ^ (Ei) . Since the sum of these level densities for the full spectrum is not 
normalized to unity, we define normalized level density as 

= 1 DH»)(£.), (5) 

[m + n)" 

such that £^X> (m|n) (£i) = 1. In the following, we shall use this V { ~ m \ n) {E i ) to study 
some spectral properties of the supersymmetric Polychronakos spin chain. 

Let us investigate at first whether, in analogy with the case of many other quantum 
integrable spin chains with long-range interaction, the energy level density of su(m\n) 
supersymmetric Polychronakos spin chain can be described by the Gaussian distribution. 
The form of normalized Gaussian distribution is given by 



G(E) = -±=- exp 

V Z7TCX 



(g ~ 

2a 2 



(6) 



where /i and o denote the mean value and standard deviation respectively. For the case 
of supersymmetric Polychronakos spin chain, these parameters are evidently related to 
the Hamiltonian 7i( m l«) in Eq. (TjQ) as 

(m + n) [m + n) 

In order to compare T)( m \ n '(Ei) with G(Ei), we have to compute the parameters /j, and a 
in terms of m, n and N. For this purpose, let us consider some trace formulas involving 
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supersymmetric exchange operators like 



tr 


- (j pMn) )2 - = 


tr [1] = 


= s N 


tr 


p(m\n) p (m|n) 
r ij r il 


= tr 


ij 


tr 


p(m]n) p(m|n) 
r ij r kl 


= s N 





tr 



ji 



P; 



(m|n) 



tr 



s N -h. 



p(m\n) p{m\n) 
ij kj 



JV-2 



where s = m + n, t = m — n and k, I are all different indices. With the help of these 
trace formulas, \i and a in Eq. ([7]) can be expressed as 



2s 2 



N 1 

i,j=l ^ij 
(#j) 



4 4.2 N 

s 4 - r 



(7 



2s 4 



r 4 + 

ij=l 



jV 

E 

i ,j,k=l 
{i+j+k+i) 



2 2 ' 
x ij x jk 



(10) 



where These two expressions encompass different types of summations 

involving the zeros of the N-th order Hermite polynomial. It is well known that, the 
zeros of the N-th order Hermite polynomial satisfy an identity of the form [40J 



iV 



^2 — = x i 

k=l lk 
(k^i) 



Starting from this simple identity, it is possible to derive many other more complicated 
form of identities involving the zeros of the Hermite polynomial. In particular, one can 
show that [40 



N 1 

< T 2 

i,j=l V 
(i+j) 



\n{n-\) 



N 



T 4 

(<#) 



N(N-l)(2N + 5). (12) 

6b 
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By using the identity (jTTj) . we also find that 

£ -^ = §^-1)^-2). ( 13 ) 

Substituting the identities given in Eqs. ffT2l) and ffTBl to Eq. ffTUj) . we finally obtain /i and 
a in terms of m, n and N as 



/i = ^(s 2 -t)iV(iV-l) 



,2 



4s 

s 4 - t 2 2fs 2 - t 2 ) 

iV(iV- l)(2JV + 5) + j — -N(N — 1)(N — 2). (14) 



72s 4 v yv y 9s 



From the above expression it is evident that, /x and a are of the order of iV 2 and N 3 ^ 2 
respectively at N — > oo limit. Interestingly, similar type of large N behaviour of these 
parameters have also been observed for non-supersymmetric Polychronakos spin chain 
associated with the Aw-i ro °t system |38j and non-supersymmetric as well as supersym- 
metric Polychronakos spin chains associated with the BC^ root system [MIES]- 

Now we are in a position to estimate the agreement between V^^^Ei) and G(Ei) 
for the case of su{m\n) supersymmetric Polychronakos spin chain. As an example, let us 
consider the simplest case of sw(l|l) Polychronakos spin. By using Mathematica, we can 
easily express the corresponding partition function (given in Eq. ([3]) with m = n = 1) as 
a polynomial of q up to iV f» 50 and obtain T>^ 1][1 \Ei) for all values of E^. By using such 
data for the specific case of iV = 30 lattice sites, we plot V^^Ei) versus through the 
dotted curve in Fig. 1. Next, by substituting iV = 30 and m = n = 1 in Eq. flHj) , we 
get fi = 217.5 and a = 46.25. The continuous curve in Fig. 1 plots G(E) as a function 
of E for these values of \x and a. Comparison between the continuous and dotted curves 
in Fig. 1 confirms that the normalized energy level density obeys Gaussian distribution 
with a very high degree of accuracy. We also calculate the mean square error (MSE) for 
this case and find it as low as 4.74 x 10 -4 . Moreover, it is found that this MSE decreases 
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gradually with the increase of the number of lattice sites. We have compared T)^ n \Ei) 
with G(Ei) for many other values of m, n and N, and always found excellent agreement 
between them for sufficiently large number of lattice sites. 

Next, our aim is to study the distribution of spacing between consecutive energy 
levels for the case of supersymmetric Polychronakos spin chain. To this end, we define 
cumulative level spacing distribution as 



where p(x) is the probability density of the normalized spacing x between consecutive 
(unfolded) energy levels of the Hamiltonian. Unfolded energy levels are obtained by ex- 
ploiting a mapping which smooths out the 'raw ' spectrum containing local fluctuations 
jH]. To define such unfolding mapping, at first it is needed to decompose the cumula- 
tive energy level density as the sum of a fluctuating part and a continuous part. This 
continuous part of the cumulative energy level density (denoted by T)(E)) is then used to 
transform each energy E iy i = 1,2, ... ,1, into an unfolded energy Si : Ei — > Si = 7}(Ei). 
Since we have found that, the energy level density of su(m\n) Polychronakos spin chain 
follows the continuous Gaussian distribution with very good approximation, f](E) can be 
expressed through the error function as 



Finally, the function p(si) is defined as the density of normalized spacing s« = {£i+\ — 
Si)/ A, where A — (Si — Si) /(I — 1) is the mean spacing of the unfolded energy levels. 

We have already mentioned that, according to a well-known conjecture of Berry and 
Tabor, the density p(s) of normalized spacing for the case of a quantum integrable system 
should obey the Poisson's law: p(s) = e~ s |39j. However, some recent works reveal that a 
class of quantum integrable HS and Polychronakos like spin chains exhibit non-Poissonian 





(15) 
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distribution of spacing between consecutive energy levels and thus violate the Berry- Tabor 
conjecture [33-38,25]. It is also found that, for the above mentioned type of spin chains, 
P(s) obeys a certain 'square root of a logarithm' law, which can be derived analytically by 
assuming a few simple properties of the related spectra. More precisely, if any spectrum 
E min = Ei < ■ ■ ■ < Ei = E max obeys the following features: 

(i) the energy levels are equispaced, or, equivalently E i+ i — Ei = A for 1 < i < I, 

(ii) the normalized energy level density is approximately given by the Gaussian distribu- 

tion, 

(iii) two tails of the level density distribution are well spread, i.e, E max — /i, \i — E min » a, 

(iv) E min and E max are approximately symmetric with respect to the mean energy /x, 
namely \E max -\- E m ^ n 2/i| << E max E m i nj 

then P(s) is approximately given by an analytic expression of the form 



p( a ) = i-- ? =^- A /M — ). ( 16 ) 

'71 S r 



where s max denotes the maximum normalized spacing, which can be estimated with great 
accuracy through the relation 



E — E 

J -'max J -'rmn 



27TCT 
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Studying the case of su{m\n) super symmetric Polychronakos spin chain for a wide 
range of values of m, n and N, we also find that the spacing between consecutive energy 
levels does not follow the Poisson distribution. So it is natural to explore the applicability 
of expression (FTO]) for the present case. To this end, let us first check whether the energy 
spectrum of su(m\n) supersymmetric Polychronakos spin chain satisfies the four condi- 
tions that have been described in the previous paragraph. We have already found that 
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this spectrum is equally spaced with unit interval and the level density follows Gaussian 
distribution with good approximation for sufficiently large values of N. So this spectrum 
evidently satisfies the first two conditions. From Eq. f[T4"|) and the expressions of lowest 
and highest energy levels given by E min = and E max = N(N — l)/2 respectively, it 
follows that both (/x — E m i n )/a and (E max — [i)jo vary as \CR when N — > oo. Therefore, 
this spectrum conforms to the third condition. At iV — > oo limit, it is also easy to find 
that 

\Pmin -\~ E max 2/i| A(t71, Tl)(E max £/ m j n ), (18) 

where X(m,n) = fei^p ■ Note that X(m, n) becomes zero for m = n and takes a finite 
nonzero value for m ^ n. Consequently, the fourth condition is satisfied only for the case 
m — n. However it can be shown that, if one drops this forth condition, Eq. ( fTBl) still 
holds within a slightly smaller range of s [35]. Therefore, it is natural to expect that P(s) 
would follow the analytical expression P(s) ( TL6l) in the case su(m\n) Polychronakos spin 
chain for all possible values of m, n and sufficiently large values of N. 

To support the above conclusion, we study numerically the nature of P(s) with dif- 
ferent values of m, n and N. For example, we may consider the particular case of sw(l|l) 
spin chain with N = 30 lattice sites. For this case, P(s) and P(s) are drawn as the dotted 
line and the continuous line respectively in Fig. 2. From this figure it is evident that, 
the cumulative distribution of spacing matches with P(s) extremely well. The MSE for 
this case is obtained as 2.628 x 10 -6 . Next, we consider the case of su(2\l) spin chain 
with N = 30 lattice sites. The corresponding P(s) and P(s) are plotted as the dotted 
line and the continuous line respectively in Fig. 3. Again, a very good agreement is found 
between these two lines with MSE given by 4.724 x 10~ 4 . By studying such particular 
cases we conclude that, similar to the case of other quantum integrable spin chains with 
long-range interaction that have been studied so far, the cumulative distribution of spac- 
ing between consecutive energy levels of the supersymmetric Polychronakos spin chain 
follows the expression ( fTBl) with remarkable accuracy. 
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Figure 1: Plot of the degeneracies V^^Ei) versus energies Ei (dotted line) for the case 
of Polychronakos spin chain with N = 30, and its comparison with the Gaussian 

distribution (continuous line). 
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Figure 2: The dotted curve represents the cumulative distribution of spacing between 
consecutive energy levels for su(l|l) spin chain with N = 30, and the continuous curve 
represents the corresponding P(s). 
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Figure 3: The dotted curve represents the cumulative distribution of spacing between 
consecutive energy levels for su(2|l) spin chain with N = 30, and the continuous curve 
represents the corresponding P(s). 
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